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Abstract: 

We study the scattering dynamics of an n-componcnt spinor wavcfunction in a random environment on 
a two-dimensional lattice. If the particle-hole symmetry of the Hamiltonian is spontaneously broken 
the dynamics of the quantum particles becomes diffusive on large scales. The latter is described by a 
non- interacting Grassmann field, indicating a special kind of asymptotic freedom on large scales in d = 2. 
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1 Introduction 

Conventional wisdom is that a classical approach of a randomly scattered particle leads to diffusion. Dif- 
fusion in quantum systems can cither be caused by particle-particle collisions or collisions with (static) 
impurity scatterers. If the latter arc randomly distributed, however, this may lead to Anderson local- 
ization rather than to diffusion (TJ [2] . This effect is particularly strong in low-dimensional systems such 
as in two-dimensional graphene sheets. The scaling approach to generic random scattering [2] indicates 
that diffusion is suppressed by Anderson localization for dimension d < 2. On the other hand, ballistic 
motion can also be ruled out, even for a finite system with random scattering [4]. It was pointed out 
by Kaveh, however, that diffusion cannot be obtained in random-phase approximation applied to t he 
disordered system [3]. 

Inspired by the recent observation of metallic behavior (i.e. diffusive or even ballistic transport) in 
disordered two-dimensional systems (graphene) [5J [S] , a general discussion of a diffusive quantum particle 
is required, which takes into account a spinor structure of the wavefunction. There are two possibilities, 
ballistic transport for finite systems [71 [5] or diffusive transport for infinite systems Here we will 
focus on infinite systems and study a quantum n-component spinor particle on a two-dimensional lattice 
with particle-hole symmetry. It will be shown that on large scales the particle diffuses on the lattice with 
d = 2. This work presents a generalization of the idea that a spontaneously broken supersymmetry can 
lead to diffusion in a system with particle- hole symmetry [9l 110) . 

The motion of a quantum particle is characterized by the transition probability P r r i (ie) for an n- 
componcnt spinor particle at site r' that moves to site r with frequency ie: 

P r , T ,(ie) = \ with K r<I ,(ie) = (Tr„ [<?,,,, (ie)Gj, T (ie)]) v = (Tr„ [G r , T ,(ie)G T , , r (-ie)\) v , 

(1) 

where G(ie) = (ie+H) 1 is the one-particle Green's function of the Hamiltonian H and (...)„ is the average 
with respect to some random scatterers. Tr„(...) is the trace with respect to the n spinor components. 
The last equation in Eq. (JTJ) follows from the Hermitean Hamiltonian: = H. 

After Fourier transformation of the two-particle Green's function K T y (ie) —> k T y (t) we study the 



motion of the quantum particle with the mean-square displacement of the coordinate 

(rl) = ■ (2) 

This expression grows linearly with time t in the case of diffusion. 

2 Model 

We consider an n-component spinor wavefunction described by the Hamiltonian matrix 

H = H o + vHi , H = {h T>a -r', a ') , Hi = {h a ^5 T ^) , v = (v r S a ^6 r ^) , (3) 

where r, r' are coordinates on the two-dimensional lattice and a, a' = 1,2, ...,n refer to the n spinor 
components. v r is a random variable with an uncorrelatcd Gaussian distribution: (v r ) v = 0, (v r v r i) v = 
gS r ,r'- In the following we assume that the Hamiltonian satisfies the generalized particle-hole symmetry 
Hj — > —UH*U^ = Hj (j=0,l), which belongs to class D according to Cartan's classification scheme [TTj . 
In terms of the Green's functions, this transformation provides a sign change of the frequency: 

G{ie) -> -UG T (ie)U t = G(-ie) , (4) 

since Hj = Hj implies Hj = —UHjU^ ( T is the matrix transposition). The Green's functions G(ie) 
and the transposed Green's function G T (ie) can be expressed in a functional- integral representation of a 
free field complex (boson) field 4>l k and a Grassmann (fermion) field §\ k , respectively. This allows us to 
construct the Bose-Fermi functional integral [12] 

(f(4>)U = J f{^)e- S V[4>] (5) 

which is normalized: 

J e- s V[<j)} = 1 . (6) 

The action S is 

5 = -^(0•(i^ o +^e)0) + . g (0•i^ 1 0) 2 (e > 0) , (7) 

with respect to the boson-fermion vector field <j> = (</>J k ,(f>^ k ) (k = 1, 2, n) and with the block-diagonal 
Hermitean matrices Hj = diag(Hj, Hj). After averaging over the random variables v T we can write 

(Gr,k-,r>,l( ie ) C Z',m-,T,n( ie ))v = ~ ($\' ,k$l ,J>1< ,-m) <t> ( 8 ) 

with (...),/) — J ...e~ s T>[(j>]. The normalization can easily be seen by performing the <f> integration before 
averaging over v r . 

An integral of the form ([6]) describes a supersymmetric field theory, meaning that it is a field theory 
for bosons as well as fermions which appear with the same Green's functions [12] . However, it should be 
noticed that supcrsymmetry is sufficient for the normalized integral but not necessary [3] . In the present 
case the boson and the fermion Green's functions arc different, provided that H T ^ H . The choice of 
different Green's functions in the action ([7]) has profound consequences in comparison with the model, 
where fermions and bosons appear symmetrically with the same Green's function, because the latter is 
subject to a larger symmetry group. This will be discussed at the end of the paper. 

Using the relation in Eq. ([8]). we can write for the expression in Eq. ([T]) as 

K r , T , = (Tr„ [G r , r ,(ie)G rl , r (-ie)]) v = -<Tr„ [G r ,r> (ze)C/G^ r (ze)C/ t ] )„ 

= U ™,n U l*ri(0l>,m$l,l4'l,n'$>,n)<l> = - U mtn U* n , (<j>l, n (j>l n ,^l A )^ . (9) 

This expression will be used subsequently to study diffusion in the particle-hole symmetric system. 
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3 Summary of the subsequent calculation 



Before embarking to the detailed calculation of K TiT i, that will lead us to a simple expression for the 
functional integral on large scales |r — r'| in terms of a saddle-point approximation, a brief outlook on 
the lengthy calculation is given in this section. In a first step we identify a symmetry in terms of a 
similarity transformation with respect to the boson-fermion structure. After introducing a new field 
in the functional integral, we apply a saddle-point approximation to the latter. It turns out that the 
above mentioned symmetry creates a two-dimensional fermionic saddle-point manifold, given by a two- 
component Grassmann field (<p,(p'). For large scales this becomes a free field and provides a diffusion 
propagator. In other words, our approximation scheme allows us to prove that the Fourier components 
of K r ^ r i(ie) describe diffusion in the large distance asymptotics: 

K^ie) ~ - K . (10) 
toe + c — c q 

with finite constants b and K, which is determined by the solution of the saddle-point equation. Moreover, 
c q are the Fourier components of 

c r r , = 16Tr n \g+y t rQ2Hig- !r yQ 2 H 1 ] , 

where the Green's functions g± arc defined as 

g± = [H ±ie + 2(Q 1 ±Q 2 )H 1 }- 1 . (11) 

Qi, Q 2 are determined by saddle-point equations. 

remark: The Green's functions g± can be considered as the self-consistent Born approximation (SCBA) 
of the random Green's functions G(±ie), where Qi, Q2 are self-energies [TBI H4] . 

4 Diffusion on large scales 

In the following we derive the asymptotic form of K T r i in Eq. (jTDJ) . 

4.1 Boson- Fermion Symmetry 

Considering the block matrix 

a e 
e b 

where the elements of the matrices A, B are complex and the elements of the matrices 8, are Grass- 
mannian, we introduce the graded trace 

Trg(| ®J =TrA-TrB , 

where Tr is the conventional trace, and the graded determinant detg [9] : 

For the special matrix H = diag(H, H T ) this gives Trg(iJ) = and detg(i? + ie) = 1. Trg and detg 
have the same properties as the conventional trace and determinant. In particular, we have the relations 
detg(i)detg(S) = detg(is) and detg(i) = exp(Trg(logi)). 
Now we consider the special matrix 



S 



? w V q) WeG), (13) 
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where G is a Grassmann algebra (i.e. iptp' = —ip'tp). Hj and S anticommutc: 

hj = ( H ^ U ) - ( , JL, - ( ° \ 



Hjtp'W J-\<p'HjW J-y-tflTtHj 



3 ' 



where the second equation follows from the assumption that <p, tp' commute with Hj. This relation 
implies that for a global 5 (i.e., 5 is constant on the lattice) 

e § H 3 e § = Hj , (14) 

which can be considered the supersymmetry of the model defined in 10 because the transformation 
connects the fermionic and the bosonic sector of the theory. For the subsequent calculations it is useful 
to notice that with TrgS* = we have 

detg(e^) = exp hrgs) = 1 . (15) 



4.2 Functional integral with nonlinear field 

Defining the tensor field 

&k& = fi.4f (i. J" = 1. 2; k, k' = 1, n) , 
we rewrite the terms in Eq. (JTJ as 

(0 ■ H 4>) = Trg(H $), {4> ■ H^f = Trg^M^) . 

Then the identity 

gTrgiH^H^) + ^TrgK^tf^ - Q)(igH^ - Q)] = g- lr Trg(Q 2 ) - 2zTrg(Qff 1 $) 
with matrix field 

w \Q r iP r 
allows us to write the interaction term as a Q integral: 

exp [-jTrgfffjMil)] - y"exp [-^Trg^M^) - .g- 1 Trg[(z.gi7 1 4> - Q){giHx§ - Q)]\ V[Q] 

= J exp [- 5 - lr rrg(Q 2 ) + 2iTrg(Q#i$)] V[Q] . 

With the expression on the right-hand side we can perform the <f> integration in the functional integral 
of Eq. ([5]), since (f> appears only as a quadratic form in the exponent. Thus we remain with a functional 
integral over Q: 

J F(4)e- 5 P[0] = J G{Q)detg(H + ie + 2QH 1 )- 1 e- 9 ~ 1 ^ si ^V[Q] (16) 

because of 

J e - T ^ A ^V[cb} = dctg(i)- 1 . 

The determinant 

J = dctg(H Q +ie + 2QH 1 )- 1 (17) 

is the Jacobian for the transformation (f> — > Q in the functional integration. The function G can be 
obtained from F by directly calculating the integrals on both sides. This, however, is a complex task for 
a general F. Here we consider only one specific case which is sufficient for a diffusive mode: 

K r ,r' = J U m ^ n U* n ,((H 1 1 Qr')mn{Hi 1 Q r )n'l)Q 

with = j ...Je- 3 ~ lTrs ^V[Q}. Moreover, we have = -WHJJ such that 

K r^ = \ U^nU^iiH^Qr'UnitfH^UBJn^Q . (18) 

^ l,m,n,n' 
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4.3 Saddle-point approximation 

The saddle-point approximation of the functional integral (|16[) is given by a solution of the saddle-point 
equation <5qS" = with 

5'=.g- 1 Trg(Q2) + iogdctg( J tf + ie + 2Q J ff 1 ) . (19) 
The saddle-point is degenerate with respect to the similarity transformation 



e^Qoe s with Q Q 



Qo 
iPo 



This covers the entire saddle-point degeneracy because we assume here that there is no additional sym- 
metry of H. 

Qo consists of two terms, namely Qq = Qi + Q 2l where Q\ (Q2) commutes (anticommutes) with S: 

e § Q e~ S = Q 1 +Q 2 e- 2§ (20) 
which implies Trg(QQ) = 0. Then the saddle-point solution contributes to the action (|19p the two terms 

[/(Q 1 i7 1 ) T C/ t = -Qifli, U(Q 2 H 1 ) T U^ = Q 2 H 1 , (21) 
where the first (second) term preserves (breaks) the symmetry of the Jacobian. These properties imply 

UgTtf = - g+ (22) 

which is consistent with Eq. (jj). 

Inserting the expression (f2"0")) into the functional integral of Eq. (|T5| results in 

G(e 5 Q e-^)dctg \h + ie + 2QA + 2Q 2 H ie 2§ ] * V[Q'] . (23) 

This indicates that Q 2 is the order parameter for spontaneous symmetry breaking. Thus we have reduced 
the integration to the nonlinear field Q' = Q 2 exp(2S), while Qo is determined by the saddle-point 
condition. 

Now we use the identity e 2S = 2(1 — S)^ 1 — 1 and define 7+ = ^g±Q 2 Hi with the help of the Green's 
functions g± in Eq. (|TT|) to obtain for the inverse Jacobian (cf. Appendix [Aj 

j-i T-i A .(+> > <, j dct(-[ff -te + 2(Q 1 -Q 2 )ffi]) - 

J — J Met l-l- 7+9595 -997.99 + 7+957-99 with J= . 24) 

det(H +ie + 2Q Hi) 

Using the identity det(A) = exp{Tr[log(A)]}, we eventually have 

J = Jexp {-Tr [log (1 + 7+9999' - 937- V 3 ' +1+^1-^')]} • (25) 

4.4 Large-scale properties 

The spatial diagonal elements of 7+9793' — 937-95' + 7+997-95' can be written as 

(7+9999' - 957-95' + 7+9?7_95')r,r = (7+ - 7- + 7+7-)r,r95r93r + ^2 7+,r,r'7-,r',r(¥V - (Pr)<f'r i ( 2 6) 

r' 

where the first part is proportional to e: 

(7+ - 7- + 7+7-)r,r = -8ieg+g-Q 2 H 1 . (27) 
The second term can also be expressed as 

^ 7+,r,r'7-,r',r(9V ~ <£r)</>r = rf r,r' VrVr' ( 28 ) 



with 

dr,r' = #r,r' ~Y] C r " ,r> ~ Cr, r ' with Cry = Tr„ [7+, r ',r7-,r,r'] • (29) 
r" 

It should be noticed in Eq. (|28|) that the spatial diagonal elements 7±, r ,r do not contribute. Moreover, 
("f+iptp' — tpj—ip' + 7+V7- v')r.r' (r' 7^ r) has at least one spatial off-diagonal factor 7±,r,r' in each term. 
Therefore, all matrix elements (j+ipip' ~ <^7_(ys' + 7 + </?7_<y9') r . r < have at least one factor 7±. r , r ' with r' ^ r, 
except for the diagonal term in (|27p which is proportional to e. 

In the next step we analyze terms that depend on the off-diagonal elements 7±.r,r' (r' r). Under 
a change of the length scale r — > Ar on the two-dimensional lattice these off-diagonal terms scale as (cf. 
Appendix [Bj 

7±jr , r , —¥ A- 2 7 ±, r , r , (r' ^ r) . (30) 

e is an arbitrarily small parameter which should be sent to zero. This allows us to replace e — > A _2 e here. 
Moreover, products of n matrices are of order A -2 ™ because 7±.r,r' decays exponentially in space due to 
the nonzero symmetry breaking term Q 2 . Therefore, the intermediate r summations do not contribute a 
factor A. Finally, the trace scales as Tr — > A 2 Tr, and we obtain from Eq. ([25]) for the scaled Jacobian 

J -)■ Ja = Ja cxp {-A 2 Tr [log (l + A~ 2 (7+ W ' - <frf_<p' + y + tpy_(p'))] } . 

Thus the large-scale limit A ~ oo reads 

Ja ~ JAexp{-Tr (7+w' - f~f-<p' + 1+<P1- 1 P 1 )} > ( 31 ) 
which is a quadratic form of tp, <p' in the exponent (i.e., (<p, ip') is a free field). This reads with Eqs. (|27|) . 

(HI 



Ja ~ Ja exp 



- (iebS r y + d r y) 



J A exp -^2n r l,ip r tp' rl , (32) 



where b = 8Tr n [(g + g^Q 2 Hi) rr }. After Fourier transformation r — > q we obtain 

rf q = c - c q and K q = - — — ^ — . (33) 

toe + c - c q 

Returning to the functional integral in Eq. (|18|) we now have an integration over ip, ip' with 

e r = -2Q 2 iVr, e r = 2t/tQ 2 ^ 

such that 

#r,r< ~ ^^^m.n(ffr 1 Q2^) mn ^t/ i V(t/ t J ffr 1 02)„'i(¥'r^) 

= ^Tr n (UU T H- 1 Q 2 )TT n (U*U^H- 1 Q 2 )(ip r ip' r ,) 

with {(p r <p' r i) = — /tr' r/det(/c). Using the Fourier components in Eq. (|33p . the Fourier transformation of 
K r r > reads 

K 

ibe + cq — c q 

where 

K= -^--TY n (UU T H- 1 Q 2 )Ti- n (U*U^H- 1 Q 2 ) . 
det(/c)<r 

This concludes our calculation of the large-scale properties of K v v i . 



K « ~ ITT—, . (34) 
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4.5 Alternative approach: Nonlinear sigma model 



Returning to the expression in Eq. (f2"3"f , we can expand the logarithm of the Jacobian in powers of Q2 
up to second order. This approximation is referred to as the nonlinear sigma model approach which is 
believed to provide a good description of the transport properties of disordered systems Q]3 . For our 
model we derive the nonlinear sigma model for the action 



S' = lop 



detg (h + ie + 2Q\Hi + 2Q 2 H ie 2S 

= log [detg (h + ie + 2(Qi + Q 2 )Hi + 2Q 2 H 1 (e 2S - 1) 
where e 2§ - 1 = 2(5 + S 2 ). With 



Go 



g+ 

-Wg-U 



9+ 
9 T + 



we can expand the action up to second order in Q2 as S' ~ So + S" with 



S" = 4Tri 



S 2 



8Trg 



8Trg 



HxS' 



(35) 



G0Q2H1 can be approximated by a gradient operator. This gives the standard form of the nonlinear 
sigma model for the last two terms, whereas the first term contributes to the symmetry-breaking term 
which is proportional to ie. Moreover, a straightforward calculation shows that the last term vanishes 
for our model 



Trg 



(g 4 



= 



(36) 



such that only the quadratic terms in tp survive in the nonlinear sigma model. This is in agreement with 
the exponent in Eqs. (pTTj) and ((321 . 



5 Discussion 

Our derivation of K r r i in the previous section was obtained without specifying Hq, Hi of the Hamiltonian. 
This prevents us from determining Qi, Q2 here because this requires the solution of the saddle-point 
equation. It is crucial though that the symmetry breaking term Q2 represents a mass to the Green's 
functions j± such that the latter decay exponentially There is no diffusion but localization for saddle- 
point solutions with Q2 = 0, as discussed for the case of Weyl fermions in Ref. [TU] . 

We leave the determination of Qi, Q2 for specific Hamiltonians to further work and study only the 
general structure of the diffusion propagator in Eq. (|34]l . For the large-scale behavior of the latter we 
consider q ~ 

K n ~ ~ (37) 



with 



and with the diffusion coefficients 



Co 



i,3 



D, 



1 



d 2 c a 



2b dqtdqj 



q=0 



J2 r r tr 3 Tr n [g + ,o : rQ2-Hig-,r,oQ2-Hi] 

Tr «[(5+.9- ( 32^l)r.r] 
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In the isotropic case (i.e. for Dij = DSij) we have 

-y 



dq = c - c q ~ bDq-, D =- 



2b dq k 

such that the diffusion propagator reads 



I d 2 c a 



ibe + Co — c q b te + Dq z 



From the diffusion propagator we can evaluate the dynamics of the quantum walk. We apply a Fourier 
transformation from frequency e to time t and get 

and a Fourier transformation from momentum q to real space coordinates r gives 

K e ~ r2/4Dt 
# q (i) -»• K{t) = T — . 

This provides the mean-square displacement as a function of time: 



ErM*) 



2Dt . (39) 



There is a simple scaling relation between the two-particle Green's function K r $ in Eq. ([T]) and 
saddle-point expression c r0 in Eq. (|29p as 

r r 

This result can be considered as an extension of the self-consistent Born approximation to -Kr,r' • 



Example: Weyl fcrmions with random gap: n — 2, H = id x O\ + id y (X2, Hi = 03, U = <Xi, Qi = 0, 
Q2 = — 1(77/2)03, where {cj} are Pauli matrices. The saddle-point equation reads in this case |10j 

Tr 2 [{g+g-) r ,r] = g^ 1 ■ 
Inserting this in our expressions above, we obtain b — Avq/g, K /b 2 — —1/4, 

c r ,o = -4r/ 2 Tr 2 [g+,o,r.9- r,o] , r k K rfi( ie ) nE r ^ Cr .° 



4e 2 l^-k 27 re 2 

r r 

Here we have fixed the cut-off A in Eq. (|4T|) such that det(rt) = 1. The conductivity a can be calculated 
from this expression via the Kubo approach by an analytic continuation e — > iuj/2 [lOj : 



W 2 ^r2A- r ,o(-W2) = -r , (41) 



which is the well-known minimal conductivity of graphene (except for an additional degeneracy factor 
4) 0. The disorder independent conductivity reflects the wellknown fact that the conductivity can not 
distinguish between ballistic and diffusive transport of Weyl fcrmions [12] . The diffusive behavior was 
also found in recent numerical simulations by Chalker et al. [T7] and Medvedyeva et al. [IB] . 
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5.1 Broken particle-hole symmetry 

We introduce a chemical potential [i that shifts away from particle- hole symmetry point by ±// in the 
Hamiltonian 

/H + fj,a 

_ [ H — /iijq 

H T -n<7 

V H T + /!<7o 

Then we define the Green's function in analogy to G(ie) as 

G(ie) = (5 + ie)~ 1 . (43) 

The generalization of transformation matrix S in Eq. (|13[) then is 

/ (pit/ 

«_ 

ip[W o o 
V o tp' 2 W o 

which anticommutes with H: SH = —HS. This implies the symmetry transformation 

e S He S = H 

and detg(e s ) = exp(Trg^) = 1. Now we can employ the expansion of Eq. (|35[) to obtain the nonlinear 
sigma model. It turns out that the fourth-order term in S does not vanish for ji ^ 0, in contrast to the 
result in Eq. 

6 Conclusions 

We have seen that the discrete particle-hole symmetry of the Hamiltonian H — > —UH*U^ = H can lead 
to a diffusive behavior. For this result it is crucial that no additional continuous symmetry exists for 
the H. A typical realization of this case are two-dimensional Weyl-Dirac fermions with random gap [5]. 
The diffusive behavior requires a non- vanishing symmetry-breaking term Qi , which reflects spontaneous 
breaking of the symmetry in Eq. I|14p . Qi must be determined as a solution of the saddle-point equation. 
This can, depending on the specific Hamiltonian H 7 generate a complex phase diagram with metallic (i.e. 
diffusive), insulating and quantum-Hall phases (c.f. |10j). 

A central fact in Sect. 14.41 is that the saddle-point integration in Eq. (|2"3"]) is restricted to a two- 
component Grassmann field (ip, ip'). This is crucial for the derivation of the main result. The integration 
would be over a larger manifold when the underlying Hamiltonian has additional symmetries or in the 
absence of particle- hole symmetry. The latter case was briefly discussed in Sect. 15. II where we introduced 
a shift away from the particle- hole symmetry point. The integration over a larger manifold may result in 
a non-diffusive behavior. 

There is a large number of publications on the subject of disordered particle-hole symmetric Hamil- 
tonians (class D), which are based on (i) field theory (in particular, nonlinear sigma models), (ii) related 
network models and (iii) numerical simulations. A discussion with many references can be found, for 
instance, in Ref. [19) . Unfortunately, there is no simple conclusion from all the publications because 
the details of the results depend on the specific form of the Hamiltonians or the network models, the 
distribution of disorder as well as on the approximations used in analytic treatments. Moreover, the 
mapping from network models onto Hamiltonian models is only understood on an approximative level 

[2812]. 

The approach discussed in this paper, which was originally proposed in Ref. [9], offers an alternative 
to the nonlinear sigma model used in Ref. |22j . The main difference between the two approaches is 
that the former is not supersymmetric, in contrast to the latter. The reason is that we started from 
the asymmetric two-particle (Bose-Fermi) Hamiltonian H = diag(H, H T ) in the construction of the 
functional integral in Eq. ([7]), whereas Bocquct et al. used the symmetric two-particle (Bose-Fermi) 





(44) 
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Hamiltonian H = diag(H, H). This difference lias several consequences for the effective field theory of the 
average Green's functions. First, the saddle-point manifold defined in Eq. (f20|) is different from the ortho- 
symplcctic Lee group OSp(2n\2n)/GL(n\n) which generates the manifold of the symmetric approach |22j . 
Second, the massless mode is only the two-component Grassmann field (05, (p 1 ) in the asymmetric approach, 
whereas it consists of Grassmann and Goldstonc (bosonic) components in the symmetric approach. Thus 
the saddle-point integration is more complex in the latter. It was treated within a renormalization-group 
approach, which provides an ideal metallic fixed point with infinite conductivity, in contrast to our finite 
conductivity in Eq. (|41[) . Besides its technical simplicity, the asymmetric approach provides a metal- 
insulator phase diagram [10] . which agrees qualitatively with the numerically determined phase diagram 
of Chalker et al. Q7]. 
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A Jacobian 

The inverse Jacobian in Eq. ([23"]! reads 

J' 1 = detg (h +ie + 2Q 1 H 1 + 2Q 2 H x e 2 ^ = dctg (6 + it + 2Q 1 H 1 - 2Q 2 H X + 4Q 2 #i(l - 5)" 1 ) • 

(45) 

After pulling out the factor (1 — 5) 1 we get 

J' 1 = detg(l - 5)" 1 detg (h +ie + 2Q H X - [H + ie + 2{Qx - Q 2 )Hx]S) • 
The (anti) commutation relation of Q\ (Q 2 ) implies 

tup, = -(-lYQjU, u^Q 3 = -li-iyp^ u = 1,2) 

and yields 

J = dot 1(1 — 0505 ) = 7p~~ 

1 y YY n det(Fj + ie + 2iP Q H?) 
xdet(l-[H + ie + 2(Q 1 -Q 2 )H 1 }p[H -ie + 2{Q 1 -Q 2 )H 1 }- 1 (H -ie + 2Q H 1 )^(H + ie + 2Q H 1 )- 1 ) . 
In the second factor, Pj can be expressed by Qj such that 

det(Hj + ie + 2iP H?) = det(U(H£ + ie + 2iP H 1 T )U^) = det(- [H -ie + 2{Qx - Q 2 )Hi\) . 

With the identities 

[H -ie + 2{Q X - Q 2 )Hx]~ X {Ho -ie + 2{Q 1 + Q 2 )H{) = l + [H -ie + 2(Qx - Q 2 )H x ]- l ±Q 2 Hx 

=:l+Ag_Q 2 Hx 

and 

[Ho + ie + 2(Qi + Q 2 )Hx}- x (H - ie + 2{Q 1 + Q 2 )Hx) =l-[H +ie + 2(Qi + Q 2 )Hx]~ l ±Q 2 Hx 

=: 1 - 4.g+Q 2 i/ 1 
and with det(l — 0505') _1 = det(l 4- ipip') we get eventually 

t-i det(H +ie + 2Q H 1 ) 

J = dct { -[H - l e + 2(Q 1 -Q 2 )H 1 }) det[1{1 + w)] 
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x dct(l - {1 - 4[ff + it + 2(Qi + g 2 ) J ffi]- 1 g 2 i?i}^{l + 4[tf - «e + 2(Qi - Q 2 )H 1 ]- 1 Q 2 H 1 }i P ') . (46) 
Moreover, we have 

det(l - {1 - 4[H + + 2(Q 1 + Q 2 )H 1 }- 1 Q 2 H 1 }( P {1 + 4[F - ie + 2{Q X - Q 2 )H 1 ]- 1 Q 2 H 1 }^) 
= det(l - W ' + 4[F + ie + 2(Qx + Q 2 )H l ]- 1 Q 2 Hi<pipf - 4^[£T - *e + 2(Qi - Q 2 )H 1 ]- 1 Q 2 H 1 <p' 
+16[H + ie + 2(Qi + Q a )#i] _1 Q2ffi¥>[flo - + 2(Qi - Q 2 )H 1 }- 1 Q 2 H 1 p') 
= det(l - <py>' + 4g + Q 2 Hnpip' - 4<pg-Q 2 Hi(p' + l6g+Q 2 Hiipg-Q 2 H 1 ip') . 
Thus, we get for the expression in Eq. (|46|) 



i det(i? + ie + 2Q Hi) , 

J' 1 = r„ ■ T 77^ n 1 rr n det 1 + ±9+Q*Kw<P - *P9- QzHkp' + 16 g+ Q 2 g iyg _ Q 2 ffx y/ 

det(- [i/ - «e + 2(Qi - C^J-HiJ) 

B Scaling transformation 

The Green's function of the saddle-point approximation in Eq. reads in Fourier representation 

rA e lqrcosa da 



o + m + 



-grf? , (47) 



where m is an effective mass that is created by the saddle-point matrices Q\ ± Q 2 . Rescaling r — > Ar 
then gives 

3Ar = / ^V- r-5— «<fc = A" 2 / -^j^pdp ~ A"V (48) 



/ ie + m + q 2 J ie + m + p 2 /A 2 

if m <~ 1, since the integral is dominated by small p and does not depend on the cut-off A A 
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